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1. The effective theory and the chromo-magnetic operator 

We consider the classical HQET Lagrangian [[l], ||, |||] of a heavy fermion of mass 1 m, whose 
spinor we indicate with \j/ h . Keeping a four component notation with P+y^h = V^h we tnus nave 

jSf = JSf stat + j£fW + 0(l//n 2 ), (1.1) 

= W h Do¥ b , = "^(^kin + ^spin) = ^^("^ " 2^«)Vh , (1-2) 

^kin = ¥h^Vh ) ^spin = YhTAl a klYh = Wh^^Yh ■ ( L3 ) 

where Z? 2 = DkDk, Gu = 5 [%, Yi] and F^ v is the QCD field strength tensor. The spin-flavor symme- 
try of the static Lagrangian stat is broken at the 0(1 /m) by the kinetic and the chromo-magnetic 
operators. At this order only the latter is responsible for the spin interaction. In particular the 
quadratic mass splitting between the ground state pseudoscalar (PS) and vector (V) heavy-light 
mesons assumes the form 

Am 2 = My — M% s = 4A 2 + 0(A^ CD /m) . (1 .4) 

The parameter A 2 is directly related to ^ sp i n and encodes, at order 1/ra, the information upon the 
deviations from the static limit, where My = Mps, stemming from the spin-dependent interactions 



inside the heavy-light mesons. The splitting (1.4) can be rewritten in two equivalent ways 



Am 2 = 4C mag (M/A)A 2 RGI + 0(f) = 2^±^C spin (M/A)A 2 RGI + 0(f) , A = A m . (1.5) 

The coefficients C mag and C sp i n perform the matching between HQET and QCD, and are expressed 
as functions of the RGI heavy quark mass M, defined as in [Q]. They are computable in continuum 
perturbation theory, and a three-loop result is presented in Sect. f|, where a motivation for preferring 
the second form in ( |l.5[ ) is provided. The RGI parameter Xf Gl is given by 

A 2 RGI = i(B|^|B)/(B|B) , = hm o [2^- 2 (M)]^ /2fo0 <i„(M) , d-6) 

with 7o = 3/(8tt 2 ), Z7 = (11 - |M)/(167r 2 ), (1.7) 

and the zero-momentum static-light meson state |B). The operator ^f pin (Ai) is related to the bare 
operator <^ S pin by a multiplicative renormalization factor in (ju) depending on the adopted scheme 
S and a renormalization scale jJ,, whereas Z^(go) = ^^/^spta depends on the bare coupling 
only. The relation between the two renormalization factors reads 

4in(M)/Zs R p?n = = US ^) , (1-8) 

where 

U s (n) = [2b Q f s ^W^exp I J* sW dg (1.9) 

is the solution of the renormalization group equation in terms of the anomalous dimension 7 s and 
the [5 -function in the S scheme with their leading order coupling expansion coefficients Here 
<I> stands for any matrix element of ^pm, e.g. %i. 

'The details upon the heavy quark mass definition are irrelevant for the present discussion. 
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2. Non-perturbative renormalization 

We follow the general strategy of and formulate a renormalization condition for ^ S pin 
in a finite volume, which enables us to non-perturbatively compute the renormalization factor 
Z^? 1 . As we are interested in accurate simulations as well as perturbative computations we choose 
Schrodinger functional (SF) boundary conditions; see for a recent review. They induce a non- 
trivial background field, F^ v , at tree-level. This ensures a good signal in MC simulations at weak 
coupling. Further, it means that a 1-loop computation is sufficient to know the renormalization 
factor up to and including 0(g ( 2 ). Since ^ sp i n does not contain any light fermion fields, we are able 
to avoid these altogether in the definition of the correlation functions. It follows that for iVf = we 
end up with a pure gauge theory definition (with no relativistic valence quarks) and the observables 
are 0(a)-improved, once the action is. 

In a discretized box of volume L 4 we adopt Dirichlet boundary conditions in the 3-direction 
and periodic boundary conditions in all others. A natural renormalization condition is then 



Z 



SF L 2 (^(x + §6)^ sp i n (x)) = L 2 {S x {x+\0)ff^{x)) 
spinl ' (S^x+ms^x)) (£!(*+ §6)5! (*)) 



, x 3 =L/2. (2.1) 

£0=0 



The spin operator Si (x) = 1+a g Vh G ^o ( x ~ a ®) Vh ( x ~ a ®) * s introduced in order to have a non- 
vanishing trace in spin space. It is a (local) Noether charge and does not need to be renormalized. 
Wo is the same temporal parallel transporter appearing in the discretized static action [||], and 
5mw is an additive mass renormalization term, whose knowledge is not needed in the following; it 



cancels out in the ratios of eq. (2.1). 



After integrating the static quark fields out and exploiting the properties of the static propagator 
[J7|, |6|], we use the equivalence of all coordinates in Euclidean space to switch to the usual SF 
boundary conditions, corresponding to "point A" in [Q], and obtain 

7SF , r 2 <Tr(^ 3 M£i(*))> _ L 2 (Tr(^ 3 «£i(*))> _7il + V3 , , rVt 



Z sp in (L) - 



(Tr(^ 3 (x))) (Tr(^ 3 «)) 



+ 0((a/L) 4 ), (2.2) 
gQ =o o 2 — v 3 



with xq = Lq/2, Ei = iFoi(x), and Dirichlet boundary conditions in time. Here Poi(x) stands for 
the clover leaf discretization of the field strength tensor [j9p. 

Having specified the lattice setup and the renormalization condition, we introduce the step 
scaling function a sp i n (M) via 

^(/i) = a spin (f 2 (l/ i u)X n (2M) . (2.3) 
It is obtained as the continuum limit 

'"^spin 



Zgn(2L) 

a/L-fO Z^- (L) 

1 spin V. / 



Cspin(«) = lim E spin (w 5 a/L) of I spin (w,a/L) 

spin 

2 



(2-4) 

g 2 (L)=u,m=0 



where g (L) is the SF coupling and the condition m = of vanishing light quark masses plays a 
role only in case that the computation is extended to Nf > 0. We performed pure gauge theory 
simulations to determine £ sp i n for different couplings u and resolutions a/L. The continuum limit 
results (see Figure |l|) allow us to reconstruct the non-perturbative scale dependence of the SF 
renormalized chromo-magnetic operator. 
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Figure 1: Left: examples of continuum limit extrapolations of E sp j n (cf. eqs. (|2.4|)) for couplings u = 
1.243,2.77 and3.48. Filled symbols indicate that F^y was defined as F^y with the link variables replaced by 
HYP2 |S] links. Right: scale dependence of ^ S pin in the SF scheme with its associated A parameter [pL 



By applying eqs. (1.7, L8|) at weak coupling g 2 (/J.) with the two-loop anomalous dimension 
in the SF scheme [[R], [TTJ] , 

r SF (^) = -^ 2 ( 7o + r f F ^ 2 + ...), yf F = -0.00236 -0.00352Af f + 0.00023Af f 2 , (2.5) 

we are able to non-perturbatively connect the low energy regime with the RGI, and arrive at 



RGI 

spin 



0.992(29), at/i = l/2L n 



2L n 



1.436r Iff 



(2-6) 



The latter has to be combined with values of Z sp j n (2L max ), depending on the bare coupling and 
lattice action, to form 

Zfcj = Zf* D (L)<t>™/^ m (l/L) (2.7) 
for the respective action. The numerical values are well represented by 



Zg tt (2L nMK ) = 2.55 + 0.16(0 - 6) - 0.40(j3 - 6) 2 , 6.0 < j3 < 6.5 , 
for the HYP1 [Hi action with an error of about 1%. For the other actions see [fT~ 



(2.8) 



3. Three-loop matching between HQET and QCD 

As pointed out in Sect, [j] the perturbative matching between HQET and QCD plays a very 
important role in a precise determination of the mass splitting. Our three-loop computation [|l3|] of 
the matching coefficient and the anomalous dimension of the chromo-magnetic operator allow us 
to give a reliable final result and estimate its uncertainty. 

The coefficient of the chromo-magnetic term needs to be determined by matching to QCD. In 
perturbation theory we consider the scattering amplitude of an on-shell heavy quark in an external 
chromo-magnetic field, expanded in the momentum transfer q up to the linear term. Denoting it 
schematically by and indicating only the presently relevant dependences, we have the (tradi- 
tional) matching condition (with =0^(aO = U MS (n) .s/^ as in eq. (|L9|)) 



1 



RGI 
hqet 



sRGI 
spin 



a) 



(3.1) 
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By working in the MS scheme and with the background field method [14], we arrive at the 3-loop 
result for the matching coefficient 



C C m(m Q ) = l+0.6897a M s(m Q ) + (2.2182 - 0.mSN f )a^(m Q ) 

+(1 1 .0763 - 1 .7495^' + 0.05 13^)0^ (m Q ) + 0(<4s 

while for the anomalous dimension of 0^ n , which enters U MS (n), we extract 



(3.2) 



7^(0^5) = 0.47750^ + (0.4306 - 0.0549^Vf)c^ 



MS 

+(0.8823 - 0.1472A& - 0.0007/^)0^ + 0(a^ 7 



(3.3) 



Here, formulae are given for the case where the heavy quarks are quenched also in QCD. Their 
loop effects are very small [13]. The conversion function C mag of Sect. [j]is obtained by changing 



the renormalization scheme in the effective theory such as to include the finite renormalization C cm , 
while C S pi n is constructed by replacing in addition the pole mass, wiq, by the RGI mass, M: 

^,cd = ^(Jf/AsgX = ^Q pm (M/A Ws )^. (3.4) 

The resulting equations 

C spin (M/A Ws ) = V&fa) = —C mag (M/A Ws ) = —U™z(m*) (3.5) 

then define the anomalous dimensions y spin , y ma §. In all these schemes the renormalization of the 
coupling remains untouched: MS. The change from the MS -mass at its own scale as argument 
of U spm to the RGI-mass as the argument of C sp i n is convenient since the RGI-masses are the 
primary quantities obtained in a non-perturbative lattice computation [ffl|. 



The second equation in (3.4) avoids the pole mass which is known to have a bad perturbative 
expansion in terms of short distance masses (or M). Thus the anomalous dimension y^ 111 is expected 
to show a better behaved perturbative series which will be reflected in C sp i n . 

For practical purposes we parametrize the conversion functions C sp i n and C mag in the Nf = 
theory, graphically represented in Figure ^[ in terms of the variable x = l/ln(M/Aj^): 

r^ p, 7PM { i + o.087x-0.02U 2 } 2-loopy «»n__ 9/ ,,o n - 

\ x iT/P«o){l+0.097* + 0.115x 2 -0.038x 3 } 3-loop y' % ~ A ] ' ( J 

These formulae guarantee at least 0.3% precision for x < 0.6. Inspection of Figure |^ shows the ex- 
pected bad perturbative behavior of C mag . We thus focus our attention on C sp i n which exhibits very 
small higher order contributions in the b-region. The difference AC sp i n (Mb/A) as 10 -2 between 
the three-loop and the two-loop determination with Mb = 6.76(9) GeV (from [|l5|]) is much smaller 
than the statistical error on the spin splitting presented in the following section. Evaluating it with 
an estimate (where the four-loop term in the very well behaved y MS is neglected) for the anomalous 
dimension y spm gives AC sp i n (Mb/A) pa 10~ 2 with respect to the three-loop estimate. We thus claim 
an about 1% relative error for C sp i n evaluated with the three-loop y spin for B-physics applications. 
For Nf = 4 the behavior of C mag and C sp j n is very similar to Fig. |2] JT3[]. 



5 



The B-meson mass splitting 



D. Guazzini 




0.8 



Figure 2: Conversion functions for Nf = 0. Dotted, dashed and solid lines use the one-, two- and three-loop 
anomalous dimension. The abscissae of the b- and c-quark [IT], [l6|, |l7[| are marked by dotted lines. 



4. First results for the spin-splitting and outlook 

As a first application we take quenched results for the bare A2 from the literature and exploit 
our results (|2.7|, ET(|). Unfortunately they exist only for j3 = 6.0, corresponding to a m 0.1 fm, 



Ref. H 

Ref. nn 



Am 2 = 0.28(6)(?)GeV 2 
Am 2 = 0.36(4) (?)GeV 2 



(00) 



Am 



Am 



0.38(7) (?)GeV 2 , (4.1) 
0.53(6) (?)GeV 2 , (4.2) 



where the numbers on the l.h.s. are taken from the corresponding references, performing a per- 



turbative renormalization. On the r.h.s. we used the b-quark mass from J15Q and the 3-loop de 



termination of C sp i n . The uncertainty marked as (?) refers to lattice artefacts and the missing 
dynamical quark determinant. The central values are now closer to the experimental mass split- 
ting, Am 2 = 0.497 GeV 2 , but at the moment the large uncertainties prevent us from concluding that 
indeed the quenched approximation can give a good estimate of this observable. 

As explained in flTo|], the same renormalization factor applies to spin-dependent potentials 



[20, 21], where so far only a perturbative renormalization was possible. 



The non-perturbative computation of Z^? 1 has demonstrated the applicability of the Schrodinger 



functional renormalization programme [22, 23] to another difficult case. Quite significant devia- 
tions from the perturbative scale evolution are present at low energies, see Figure [l[ 

With respect to a perturbative estimate, the new Zf^ has a rather big effect. Furthermore, 



thanks to the results presented in Sect. [3[ which extend [|^, 25, 26, 27, 28], we can match the 
effective theory and QCD introducing an error in practice negligible in comparison to all other 
uncertainties entering Am 2 . It now remains to compute X\ w<t with higher precision and perform 
the continuum limit. However, due to the large amount of statistics needed especially at large cou- 
plings, an extension of this method to the dynamical quarks case seems difficult. In this direction, 
other, fully non-perturbative, approaches are more promising at present fll6| , O, 29, 15]. 
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